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Abstract: Dwarf Spheroidal galaxies are amongst the best targets to search for a Dark Matter (DM) 
annihilation signal. The annihilation of WIMPs in the center of Sagittarius dwarf spheroidal (Sgr dSph) 
galaxy would produce high energy 7-rays in the final state. Observations carried out with the H.E.S.S. 
array of Imaging Atmospheric Cherenkov telescopes are presented. A careful modelling of the Dark 
Matter halo profile of Sgr dwarf was performed using latest measurements on its structural parameters. 
Constraints on the velocity-weighted cross section of Dark Matter particles are derived in the framework 
of Supersymmetric and Kaluza-Klein models. 
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Introduction 

Astrophysical and cosmological observations pro- 
vide a substantial body of evidences for the ex- 
istence of Cold Dark Matter (CDM) although its 
nature remains still unknown. It is commonly as- 
sumed that CDM is composed of yet undiscov- 
ered non-baryonic particles for which plausible 
candidates are Weakly Interacting Massive Parti- 
cles (WIMPs). In most theories, candidates for 
CDM are predicted in theories beyond the Standard 
Model of particle physics [6]. The annihilation of 
WIMPs into 7-rays may lead to detectable very 
high energy (VHE, E > 100 GeV) 7-ray fluxes 
above background via continuum emission from 
the hadronization and decay of the cascading an- 
nihilation products, predominantly from tt'^'s gen- 
erated in the quark jets. Among the best-motivated 
CDM candidates are the lightest neutralino x pro- 
vided by R-parity conserving supersymmetric ex- 
tensions of the Standard Model [12], and the light- 
est Kaluza-Klein particle (LKP) [19] in universal 
extra dimension theories which is most often the 
first KK mode of the hypercharge gauge boson. 



B'^'^\ The H.E.S.S. array of Imaging Atmospheric 
Cherenkov Telescopes (lACTs), designed for high 
sensitivity measurements in the 100 GeV - 10 TeV 
energy regime, is a suitable instrument to detect 
VHE 7-rays and investigate their possible origin. 
Dwarf Spheroidal galaxies such as Sagittarius 
or Canis Major, discovered recently in the Lo- 
cal Group, are among the most extreme DM- 
dominated environments. Indeed, measurements 
of roughly constant radial velocity dispersion of 
stars imply large mass to luminosity ratios [20]. 
The core of the Sgr dSph at 1=5.6° and b=-14° in 
galactic coordinates at a distance of about 24 kpc 
from the Sun [14]. Latest velocity dispersion mea- 
surements on M giant stars with 2MASS yields 
a mass to light ratio of about 25 [13]. The lu- 
minous density profile of Sgr dSph has two com- 
ponents [16]. The compact component, namely 
the core, is characterized by a size of about 3 pc 
FWHM, which corresponds to a point-like region 
for H.E.S.S. This is the DM annihilation region 
from which 7-ray signal may be expected. A dif- 
fuse component is well fitted by a King model with 
a characteristic size of 1.6 kpc. 
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We present in this paper the observations of the Sgr 
dSph galaxy by the H.E.S.S. array of Imaging At- 
mospheric Cherenkov Telescopes. A careful mod- 
eling of the Dark Matter halo using the latest mea- 
surements on the structural parameters of Sagittar- 
ius is presented to derive constraints on the WIMP 
velocity- weighted annihilation rate. 

Search for VHE 7-rays from observa- 
tions of Sagittarius dwarf by H.E.S.S. 

H.E.S.S. (High Energy Stereoscopic System) has 
observed the Sgr dSph in June 2006 with zenith 
angles ranging from 7° to 43° around an average 
value of 19°. A total of 11 hours of high qual- 
ity data are available for the analysis after standard 
selection cuts. After calibration of the raw shower 
images from PMT signals [1], two independent 
reconstruction techniques were combined to se- 
lect 7-ray events and reconstruct their direction 
and energy. The first one uses the Hillas moment 
method [2]. The second analysis referred here- 
after as "Model Analysis", is based on the pixel- 
per-pixel comparison of the shower image with a 
template generated by a semi-analytical shower de- 
velopment model [7, 8]. The separation between 
7 candidates and hadrons is done using a combi- 
nation of the Model goodness-of-fit parameter [8] 
and the Hillas mean scaled width and length pa- 
rameters, which results in an improved background 
rejection [3]. An additional cut on the primary in- 
teraction depth is used to improve background re- 
jection. 

The on-source signal is defined by integrating all 
the events with angular position in a circle around 
the target position with a radius of 9 cut- The tar- 
get position is chosen according to the photomet- 
ric measurements of the Sgr dSph luminous cusp 
showing that the position of the center corresponds 
to the center of the globular cluster M 54 [15]. 
The target position is thus found to be (RA = 
18''55'"59.9s, Dec = -30d28'59.9") in equato- 
rial coordinates (J2000.0). The signal coming from 
Sgr dSph is expected to come from a region of 
1.5 pc, about 30", much smaller the H.E.S.S. point 
spread function (PSF). A 6cut value of 0.14° suit- 
able for a point-like source was therefore used in 
the analysis. In case of a Navarro-Frenk- White 
(NEW) density profile [17] for which p follows r^^ 



or a cored profile [9] folded with the point spread 
function (PSF) of H.E.S.S., the integration region 
allows to retrieve a significant fraction of the ex- 
pected signal. See Table 1. 

No significant 7-ray excess is detected in the sky 
map. We thus derived the 95% confidence level 
upper limit on the observed number of 7-rays: 
]\j95%c.L. rpjjg jimit is computed knowing the 
numbers of events in the signal and background 
regions above the energy of 250 GeV using the 
Feldman & Cousins method [10] and we obtain: 
]\j95% c.L. _ gg Given the acceptance of the de- 
tector for the observations of the Sgr dSph, a 95% 
confidence level upper limit on the 7-ray flux is 
also derived: 

$^(E^ > 250 GcV) < 3.6 x lO^^^cm^^g-i 

(95% C.L.) (1) 



Predictions of 7-ray from Dark IMatter 
annihilations 

The 7-ray flux from annihilations of DM particles 
of mass niDM accumulating in a spherical DM 
halo can be expressed in the form: 
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as a product of a particle physics component with 
an astrophysics component. The particle physics 
part contains {<tv), the velocity-weighted annihila- 
tion cross section, and dN^/dEj, the differential 
7-ray spectrum summed over the whole final states 
with their corresponding branching ratios. The as- 
trophysical part corresponds to the line-of-sight- 
integrated squared density of the DM distribution 
J, averaged over the instrument solid angle integra- 
tion region for H.E.S.S. (An = 2 x lO^^sr): 

J{An) = -t'tt / PSF* / p'^{r[s])dsdn (3) 
All J An Ji.o.s 

where PSF is the point spread function of H.E.S.S. 
The mass distribution of the DM halo of Sgr dwarf 
has been described by plausible models taking into 
account the best available measurements of the Sgr 
dwarf galactic structure parameters. We have used 
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Halo type 


Parameters 


J 

(lO^^GeV^cm-S) 


Fraction of signal 
in Arj = 2 X lO^^sr 


Cusped NFW halo 


r^ = 0.2 kpc 

A = 3.3 X 1O^M0 


2.2 


93.6% 


Cored halo 


Yc = 1.5 pC 

Va = 13.4 km s"^ 


75.0 


99.9% 



Table 1: Structural parameters for a cusped NFW (rg, A) and a cored (re, Va) DM halo model, respectively. 
The values of the solid-angle-averaged l.o.s integrated squared DM distribution are reported in both cases 
for the solid angle integration region Ail = 2 x lO^^sr. 



two widely different models. The first has a NFW 
cusped profile [17] with the mass density given by: 
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with A the normalization factor and rs the scale 
radius taken from [9]. Using Eq. 3, the value of J 
obtained with this model is reported in Table 1 . We 
have also studied a core-type halo model as in [9] 
characterized by the mass density: 
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where r^ is the core radius and Va a velocity scale. 
However, we have tried to update the Va and r^ 
values which were used in [9]. By inserting in 
the Jeans equation the luminosity profile of the Sgr 
dwarf core of the form: 
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we estimated from the data of [15] a — 2.69±0.10 
and Tc = 1.5 pc. Note that the value of Tc is 
only an upper limit. The value of the central ve- 
locity dispersion of Sgr Dwarf is cr = 8.2 ± 
0.3kms~^ [21]. We have assumed that the veloc- 
ity dispersion is independent of position. The value 
of Va is then given by Va = \/a cr = 13.4 km s^^. 
The cored model gives a very large value of J, 
which is reported in Table 1 . The third column of 
Table 1 gives the amount of signal expected in the 
solid angle integration region Ai7 = 2 x lO^^sr 
Fig. 1 shows the limits in the case of a cored (green 
dashed line) and cusped NFW (red dotted line) pro- 
file using the value of J computed above. Pre- 
dictions for phenomenological MSSM (pMSSM) 



models are displayed (grey points) as well as those 
satisfying in addition the WMAP constraints on the 
CDM reHc density (blue points). The SUSY mod- 
els are calculated with DarkSUSY4.1 [11]. In the 
case of a cusped NFW profile, the H.E.S.S. ob- 
servations do not establish severe constraints on 
the velocity-weighted cross section. For a cored 
profile, due to a higher central density, stronger 
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Figure 1: Upper limits at 95% C.L. on {av) ver- 
sus the DM particle mass in the case of a cusped 
NFW (red dotted line) and a cored (green dashed 
line) DM halo profiles respectively. The pMSSM 
parameter space was explored with DarkSUSY 
4.1 [11], each point on the plot corresponding 
to a specific model (grey point). Amongst these 
models, those satisfying in addition the WMAP 
constraints on the CDM relic density are overlaid 
as blue square. The limits in case of neutralino 
dark matter from pMSSM are derived using the 
parametrisation from [5] for a higgsino-type neu- 
tralino annihilation 7 profiles. 
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constraints are derived and some pMSSM models 
can be excluded in the upper part of the pMSSM 
scanned region. 

In the case of KK dark matter, the differential 7 
spectrum is parametrized using Pythia [18] sim- 
ulations and branching ratios from [19]. Predic- 
tions for the velocity-weighted cross section of 
B^^^ dark matter particle are performed using the 
formula given in [4]. In this case, the expression 
for {av) depends analytically on the B(i) mass 



Conclusions 

The observations of Sgr dSph with H.E.S.S. reveal 
no significant 7-ray excess at the nominal target 
position. The Sagittarius dwarf DM halo profile 
has been modeled using latest measurements of its 
structure parameters. Constraints have been de- 
rived on the velocity-weighted cross section of the 
DM particle in the framework of supersymmetric 
and Kaluza-Klein models. 
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Figure 2: Upper limits at 95% C.L. on {av) ver- 
sus the DM particle mass in the B*^^) Kaluza-Klein 
scenarios for a cusped NFW (red dotted line) and 
a cored (green dashed line) DM halo profiles re- 
spectively. The blue line corresponds to Kaluza- 
Klein models [19]. Overlaid (yellow line) are the 
KK models satisfying WMAP constraints on the 
CDM relic density. 

square. Fig. 2 shows the sensitivity of H.E.S.S. in 
the case of Kaluza-Klein models where the hyper- 
charge boson B^^^ is the LKP, for a cored (green 
sohd line) and a cusped NFW (red sohd line) pro- 
file respectively using the value of J computed in 
section 3.2. With a NFW profile, no Kaluza-Klein 
models can be tested. In the case of a cored model, 
some models providing a LKP relic density com- 
patible with WMAP contraints can be excluded. 
From the sensitivity of H.E.S.S., we exclude B*^^) 
masses lying in the range 300 - 500 GeV. 
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